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I. INTRODUCTION 

In cosmology, solutions to Einstein's field equations are obtained by imposing symmetries 
[1] on space-time. Usually, spatial homogeneity is one of the reasonable assumptions (in 
an average sense) while for cosmological phenomena over galactic scale or in smaller scale, 
inhomogeneous solutions are useful. Szekeres [2] in 1975 gave a class of inhomogeneous 
solutions representing irrotational dust. The space-time represented by these solutions 
has no killing vectors and it has invariant family of spherical hypersurfaces. Hence this 
space-time is referred as quasi-spherical space-time. 

An extensive study of gravitational collapse [3-8] has been carried out of Tolman-Bondi- 
Lemaitre (TBL) spherically symmetric space-times containing irrotational dust to support 
or disprove the cosmic censorship conjecture (CCC). A general conclusion from these 
studies is that a central curvature singularity forms but its local or global visibility depends 
on the initial data. Also over the past decades, the role of pressure within a collapsing 
cloud has been studied [9-13] but the actual role of pressure in determining the end state 
of a continual collapse is not yet clear. 

On the other hand, there is very little progress in studying non-spherical collapse. 
Basically, the difficulty is the ambiguity of horizon formation in non-spherical geometries 
and the influence of gravitational radiation. Though there is hoop conjecture by Thorne 
[14] to characterize the formation of horizon but only few works [15-19] have been done to 
confirm or refute the conjecture. Recently, an extensive study of irrotational dust collapse 
has been done in quasi-spherical Szekeres space-time both for four [20] and higher [21] 
dimensions. In this paper, we have done an extensive analysis of Szekeres model with matter 
containing pressure and studied the collapsing procedure to examine the role of pressure in 
characterizing the final singularity. The paper is organized as follows: The Szekeres model 
[22] has been described in section II. The perfect fluid solution with asymptotic behaviour 
has been presented in section III. Section IV deals with the tangential stress only while 
section V deals with gravitational collapse with non-isotropic pressure. Finally, the paper 
ends with a brief discussion in section VI. 
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II. THE SZEKERES' MODEL 

The metric ansatz for the four dimensional Szekeres' space-time [22] is of the form 

ds 2 = dt 2 - e 2a dr 2 - e 2 <\dx 2 + dy 2 ) (1) 
where a and (3 are functions of all space-time variables i.e., 

a = a(t,r,x,y), (3 = (3(t, r, x, y). 

Considering both radial and transverse stresses the energy momentum tensor has the 
following structure 

T% = diag(p, -p r , -p T , -p T ) (2) 
Hence the full set of Einstein equations are 

2d/3 + (3 2 - e- 2(3 (a 2 x + a 2 y + a xx + a yy + /3 XX + /3 yy ) + e- 2a (2a'P' - 3/3' 2 - 2(3") = p (3) 

3$ 2 +2(3- e- 2a (3' 2 - er 2 \(3 xx + (3 yy ) = -p r (4) 

a + a 2 + (3 + (3 2 + a(3 + e - 2a {a' ' (3 1 - (3 /2 - (3") - e- 2fJ (a 2 y + a yy - a y (3 y + a x (3 x ) = -p T (5) 
a + a 2 + (3 + (3 2 + a(3 + e,- 2a {a' (3' - (3 12 - (3") - e- 2( \a 2 x + a xx - a x (3 x + a y (3 y ) = -p T (6) 

a y {-a x + (3 X ) + a x (3 v ~ a xy = (7) 
a(3' -[3(3' -p' =0 (8) 

- aa x + (3a x - d x - (3 X = (9) 

— aa y + (3a y — d y — (3 y = (10) 

a x (3'-(3' x =Q (11) 

a y (3' -p v =0 (12) 

where dot, dash and subscript stands for partial differentiation with respect to i, r and the 
corresponding variables respectively (e.g., (3 X = |^). 

Combining the time derivatives of equations (11) and (12) with the differentiation of (8) 
with respect to x and y separately we have the integrability condition 

(3'(3 X = = (3'p y (13) 
Hence we may have three possibilities namely, 
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(a) p ? 0, X = = V , 

(b) 0' = 0, X = = y (14) 

(c) /?' = 0, 4 £ 0, 4 ^ 

Following Szekeres' [22] the field equations are not solvable for the third case so we shall 
consider only the first two cases. 

The energy conservation equation namely T^ v ; v — gives 

p + p(a + 20) + (ap r + 20p T ) =0 (15) 
d 

— (p T e a )=p r a x e a (16) 

Q 

— (p T e a ) = Pr a y e a (17) 

i^)=p T ^) (18) 

Now for the first choice namely 0' ^ 0, X = = y we have from the field equations the 
explicit form of the metric coefficients are as follows: 

e = R(t, r) e v{ - r ' x '^ (19) 



e a = R' + R v' (20) 
where R and v satisfy the following differential equations 

2RR + R 2 + p r R 2 = f(r), (f(r) = arbitrary separation function) (21) 

and 

e- 2 \v xx + v vv ) = f(r) - 1 (22) 
Here we have assumed p r = p r {r 1 1). Equation (21) has the first integral 

R 2 = f(r) + ^ - \ J PrR 2 dR (23) 
while one of the possible solutions of equation (22) can be taken as [24] 

e~ v = A(r)(x 2 + y 2 ) + B x {r)x + B 2 {r)y + D{r) (24) 
where the arbitrary functions A(r), B\(r), B 2 (r) and D{r) are related as 



B\ + B\ - AAD = f(r) - 1 



and G(r) is an arbitrary function. 

For the choice (b) the metric coefficients are of the form 

e f3 = R(t)e u{x - v) (25) 

e a = R(t)r ] (r,x,y)+ f x(t,r) (26) 

Then as before from the field equation (4) we have similar differential equation in R and 
v as 

2RR + R 2 + p r R 2 = K, (K is a constant) (27) 

e~ 2u {v xx + v yy ) = K (28) 
with K an arbitrary constant. Also for v we choose as in case (a) 

e~ v = P(x 2 + y 2 ) + Q x x + Q 2 y + S (29) 
where P, Qi, Q2 and S are arbitrary constants restricted by the relation 

Qi + Ql- 4PS* = K 
Now to determine the function r] we have from the field equation (7) 

and then from the field equations (5) and (6) we have a possible solution 

e~ v ri = u(r)(x 2 + y 2 ) + v\{r)x + v 2 {r)y + w(r) (31) 
where u(r), v\(r), v 2 (r) and w(r) are arbitrary functions of r alone. 
Also the differential equation in /i is 

Rji + Rfi + n(R + p T R) + (p T - p r )R 2 v = K r ) (32) 

with 

h{r) = 2{uS + wP) - (v 1 Q 1 + v 2 Q 2 ) . 
Now for explicit solutions, we shall consider the following cases in the next sections: 

(i) the perfect fluid model (i.e., p r = p T ) 

(ii) the tangential stress only (i.e., p r = 0, p T ^ 0) 

(iii) the general case (i.e., p r 7^ 0, p T 7^ 0). 

In fact, cosmological solutions are obtained (in sections III and IV) for the first two cases 
respectively while for the third case collapsing behaviour has been studied and the role of 
pressure has been examined. 
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III. THE PERFECT FLUID MODEL 

In this case due to energy conservation equations (16)-(18) the isotropic pressure is 
function of t only i.e., p = p{t) (p r = p T = p). As there is no restriction on the energy 
density so p is in general a function of all the 4 variables i.e., p — p(t, r, x, y) and hence no 
equation of state is imposed. 

Now for explicit solution according to Szafron [25] and Szafron and Wainwright [26] we 
choose 

P (t) = Pc r s (33) 

(p c and s are positive constants) and we have the general solution for R as [23] 



(34) 



I cit* + c 2 t 1 - qi 



according as £ is an integer, non-integer and s = 2. Here C\ and Ci are arbitrary functions 
of r and we have chosen 



i 



— , *="f' qi = 2 (1 + ^ ri ^ ) - 



It is to be noted that to derive the above solution we have chosen f(r) = 0. Further, 
if we consider the dust model (i.e., p = 0) then the above solution simplifies to R oc t 2 / 3 
which is the form of the scale factor in the usual Friedman model. 

Now, the physical and kinematical parameters have the following expressions 



G'(r) + 3GV _ pc 
R 2 (R' + Rv 1 ) ~ t s 



(35) 



RR! + mkv' + 2RR' 
R(R' + Rv') 



(36) 



1 
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R(Rf - 2R'f) + (RG'(r) - 3R'G) 



RR 2 (R' + Rv>) 



(37) 



The above solution is for the choice (a) (see eq.(13)). For the choice (b) (i.e., eq.(14)) the 
explicit form for R is same as in equation (34) except here C\ and C2 are arbitrary constants 
and K = 0. But we note that the differential equation (30) is not solvable for the above 
explicit solutions for R. The physical and kinematical parameters have the expressions as 



2(R 2 -K) 2{jl + r]R) p c 



R 2 



fi + r/R t s 



(38) 



Rfi + 2p,R + 3RRi] 
R(p + rjR) 



(39) 



G 



o 2 = 



Rfi — R/j, 
R{p + rjKj 



-I 2 



(40) 



Asymptotic Behaviour 

We shall now discuss the asymptotic behaviour of the above solutions. The co-ordinates 
vary over the range : to < t < oo; — oo < r < oo; — oo < x, y < oo. For both the choices 
(a) and (b) as p > 0, p^O so we must have | < qi < 1. So for s = 2 for large t, 



2<31 

R^t^ 



a 2 ~ i" 2 



(41) 



Hence as t — > oo, (p, p) fall off faster compare to (9, a), while the scale factor i? grad- 
ually increases with time. So the model approaches isotropy along fluid world line as t — > oo. 



IV. MODEL WITH TANGENTIAL STRESSES ONLY 

For this model we have from the conservation equation (18) 0' = i.e. choice (b) is only 
possible here. Also from the other conservation equations namely equation (16) and (17) we 
have 



P T = A{r,t)e~ 



(42) 



where A(r,t) is an arbitrary function of r and t. But for the consistency of the differential 
eq.(30) p T (as stated earlier) must be a function of r and t and hence a should be independent 
of x and y. As a consequence, in the solution (26) for a, r\ must be independent of x and y. 
Thus for the solution (32) for r/ we should have 

u(r) = Pr) {r), vi(r) = QiTj (r), v 2 (r) = Q 2 r)o( r ) and w(r) = Srj (r) 
Hence we have rj — r]o(r), an arbitrary function of r alone and h(r) = —Krj (r). 

Now, the differential equation for R has the simple form 



R 2 = oi + ^ , (ai, a 2 are constants) 



(43) 



which has a parametric solution of the form 



^=2ra( 1 -cos0) 
tc~t= 2( „° i 2 )3/ 2 (4> ~ sin<P) 



> for oi < (0 < < 27r) 



(44) 
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R=^(coshcf>-l) 



tr-t 



2 a 



(sinh0 



for ai > (0 > 0) 



and 



i?=(^) 1/3 (tc-t) 2/3 



for di = 



Here i c is an integration constant that corresponds to the time of arrival of each shell to 
the central singularity. 

Choosing the power law solution (i.e. a\ = 0) for R (i.e., if R oc T 2 / 3 then equation (27) 
implies that if = 0) and assuming p T — p 0T /T 2 , (T = t c — t) (i.e., a function of T alone), it 
is possible to have a solution for p (from eq.(32)) as 



' d{r)T^ + C 2 (r)T n \ p QT < 1/4 

p(r,t) = I Ci(r)T 1 / 6 cos(A: InT) + C 2 (r)T 1 / 6 sm(k InT), p 0T > 1/4 (45) 
_ C 1 (r)T 1 /6 + C 2 (r)T 1 /6 lnTi ^ = 1/4 

where Ci and C2 are arbitrary functions of r, n\ = g + — 4p 0T , n2 = g — ^^/l — 4p 0T , 

1/3 

A: = \\/^p 0T — 1 and i?o = (^f 2 ) (Here we have set 77 = 0, which has no effect on the 
metric) . 

Further, the physical and kinematical parameters have the expressions 



2(R 2 - K) 2(/i + t]R) 



R 2 



p + rjR 



2p T 



(46) 



Rp + 2pR + 3r]RR 
R(p + r)R) 



(47) 



o 2 = 



pR — pR 
R{p + r)R) 



(48) 



It is to be noted that the solution for R does not depend on p T so R has same expression 
for dust model. But for the solution of p we have only the power law form T 2//3 (or 
T -1 / 3 ) when R has Friedmann like behaviour (i.e. R ~ T 2 / 3 ). The difference comes in the 
matter density. For dust model p is a function of all the four co-ordinate variables while 
in the presence of tangential stress p is a function of t and r only. Finally, the asymptotic 
behaviour for both the model will be very similar. 



V. ROLE OF PRESSURE IN GRAVITATIONAL COLLAPSE 

In the general case when both radial and tangential pressures are non-zero and distinct 
then from the Einstein equations they can be obtained in compact form as 

Pr = (49) 
P T =Pr + ^T 
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where F(r, t) = Re 3l, {R 2 - f(r)) and ( = e^. 

Since p r is regular initially at the centre and blows up at the singularity so we can choose 
p r to be of the form: 



(») 

where g(r) is an arbitrary function such that g(r) ~ r" near r = to make initial matter 
density non-zero at the centre r = and n is any constant. As a consequence, the expressions 
for matter density and tangential stress become 



H' + ZHv' 
R 2 (R' + Rv') 



(51) 



P T = 



R n 



nR! 
2(R' + Rv') 



+ 



9'(r) 
2R n ~ 1 (R' + Rv') 



(n^3) 



(52) 



where H(R,t) = C{r) - i? 3 "" and C{r) is an arbitrary integration function. Also 

the radial velocity of collapsing shells at a distance r from the centre is given by 



R 2 = fir) + (53 ) 

Now if we choose R = r initially then at the beginning of the collapse the density and the 
tangential stress have the initial values 

... , d + 3cm' 

Pi (r, x, y) = Pi (r, t u x, y) = + (54) 



p =p T (t = U 



9(r) 



2(1 + rv') 



g'(r) 

2r n - 1 (l + ru') 



(55) 



where c(r) = H(r,U) = C(r) - f^- r 3 "" 



Here it is to be noted that for regular initial data C(r) and g(r) to be C°° functions and 
hence we have the following series expansions 



g(r) = E7= gjr n+j 

Pi(r)=ET=oPjr j (56) 

EOO -j 
j=oPj rJ 

where v_ x > — 1. 

In these scries expansions the coefficients g^s and Cj's are constants while Pj's, Vj's and 
Pj's are functions of x and y. These coefficients are related among themselves through the 
relations (54) and (55) as follows: 
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P0=90, Pi =fll{l+ 2(1+^) }' P2 = g2 {l + XT+^j} - 2(l + ,,_> ■ 

4+3!/ , 5+3i/ , 
p = 3c , pi = -^p^- Ci, P 2 = C 2 - (1+ ^). , 



(57) 



or 

Po=.9o + ^, Pi=fli(l-^)+S' P 2 = .9 2 (l-f)+ M ^p 1 .9i + if, 

according as v_ x > — 1 or v_ x = — 1 and a — d — i = 0, 1, 2, .... 

The hypersurface t — t s (r) describing the shell focusing singularity is characterized by 



(58) 



R(t s (r),r)=0 (59) 

As the differential equation in R (i.e., eq.(53)) is not solvable so we shall consider only 
the marginally bound case (i.e., f(r) = 0). Hence in this case, the singularity hypersurface 
can be written in explicit form as 

where 2-F1 is the usual hypergeometric function with b = 2 ( 3 3 _„) ■ 



A. Formation of Trapped Surfaces 



The event horizon of observers at infinity plays an important role in the nature of the 
singularity. As formation of event horizon depends greatly on the computation of null 
geodesies whose computation are almost impracticable for the present space-time geometry, 
so we consider closely related concept of a trapped surface (namely a compact space-like 
2-surface whose normals on both sides are future pointing converging null geodesic families). 
In fact, if the 2-surface S r ,t {r ^constant, t =constant) is a trapped surface then it and its 
entire future development lie behind the event horizon provided the density falls off fast 
enough at infinity. So mathematically, if if M denotes the tangent vector field to the null 
geodesies which is normal to S r j then we have 

Ky, K» = 0, K» v K v = . 

Now the convergence or divergence of the null geodesies is characterized by the sign of the 
invariant evaluated on the surface S r ,t = (in fact, < indicates convergence 

while K^. > stands for divergence). It can be shown that the inward geodesies converges 
initially and throughout the collapsing process while the outward geodesies diverges initially 
but becomes convergent after a time t a h{r) (time of formation of apparent horizon) given 

by 

R(t ah (r),r) = -,/T+W) 
Now using equations (23) and (50) we have 
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g{r)R 3 - n {t ah (r), r) - (n - 3)R(t ah {r),r) + (n - 3)C(r) = (61) 
or using the explicit solution for R (from eq.(53)) 

2r 3/2 1 ff(r)r 3 -" 1 2fl 3 / 2 (t a ,,r) 1 ff^-r^^) 

tahir) - U = WW) 2Fl [ 2 A C(r)(3-n) ] - 3^ ^ [ 2 ' *' ^ C(r)(3-n) 1 

(62) 

From the equations (60) and (62) we see that the shell focusing singularity that appears 
at r > is in the future of the apparent horizon. 

As we are interested in central shell focusing singularity (at r = 0), so its time of occurrence 
is given by 

to = t s (r) 
lira r — > 

(63) 



= U + ^2F 1 [^b,b+l,z], (z=c^)) 



where in evaluating the limit we have used the series form of g(r) and C(r) (from cq. 
(56)). Now if we restrict n < 3 then we have a comparative expression between t a h{r) and 
to as 



t a h{r) - 1 



,« |i, M + 1, „ + ffi-?$ qr" » + 1, » + a, .] 



/~i3-n i 

+ °^ - (3-n)(9-2n) ^[^M+LCg- .] r»-»" + , (n < 3) (64) 

Note that here t is the time of formation of singularity at r = while t a h{r) is the 
epoch at which a trapped surface is formed at a distance r. Thus if trapped surface 
is formed at a later instant than to then it is possible that any light signal from the 
singularity can reach an observer. As the geometry of the present model is a class 
of spherical space-time having different centres, so the condition for formation of NS 
(or BH) will be same as TBL model. Therefore, t a h{r) > to is the necessary condi- 
tion for formation of naked singularity, while to form black hole the sufficient condition 
is t a h{r) < t - It should be mentioned that this criterion for naked singularity is purely local. 

Due to complicated form of equation (64) it is very difficult to make a comparative study 
between t a h and to- Hence for simplicity we choose n = 3/2. Then the difference between 
t a h and t has the form 

2 (c gi -C\g - giVCo~\/Co - ho) 

tah(r)-t = ^ = V / r + Q(r 2 ) (65) 

ZgoVCoJ Co - 3.90 [yCo + J Co - gfiro) 

Hence in the present problem it is possible to have local naked singularity or a black hole 
form under the conditions shown in the following table (see Table I): 
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Fig.l Fig.2 



Figs. 1 and 2 show variation of t a h — to of eq.(65) for the variation of fco(= go/Co) and 
ki(— gi/Ci). Fig.l corresponds to Ci > while Fig.2 corresponds to Ci < 0. 



TABLE-I 



Choice of the parameters 


Naked Singularity 


Black hole 


(i) 5 i > 0, d < 


Always possible 


Not possible 


(ii) 9l < 0, a > 


Not possible 


Always possible 


(iii) 5i > 0, G x > 










t>i( i +V l -i fc °) 


^i+yi-ffco) 


(iV) 5i < 0, Gx < 




) 


l£l>§ 




|^|<|(l + v /l_|fc 





Here we note that for initial density gradient to be negative at the centre (i.e., px < 0) 
we must have (Ci — -^-) < (for u_ 1 > —1). In the first case (i.e., 51 > 0,Ci < 0) 
we have negative definite px and there is always naked singularity as in the dust model. 
Similarly in the second case (i.e., 51 < 0, C\ > 0), pi is positive definite and we always 
get black hole same as dust model. For the third and fourth cases (when 51 and Gx 
are of same sign) both naked singularity (NS) and black hole (BH) are possible for the 
restrictions given in the table I. When both 51 and Cx are positive (third case) or negative 
(fourth case) then for formation of NS px is negative but for BH case as there is no 
lower limit (or upper limit) of j^h (or |^-|) so px > or p\ < are possible. Further for 
gx = Cx = we have px = then we have similar behaviour for the parameters (52,(^2). 
A diagrammatic representation of t a h — to for variation of fco(= go/Co) and fci(= gx/Cx) 
has been shown in figures 1 and 2 for positive and negative Cx respectively. In both the 
figures the vertical positive region corresponds to NS while the negative region stands for 
BH solution. Finally, we see that if the initial density or pressure has opposite behaviour 
(i.e., one increases and other decreases and vise- versa) near the centre r = then we have 
similar character of the singularity as in dust model i.e., pressure has no significant effect 
on the singularity formation. On the other hand, if initial density and pressure increase 
or decrease simultaneously near the centre then even for negative density gradient at the 
centre it is possible to have a BH formation at r = 0, which is a distinct result in com- 
pare to dust model. Therefore, we may conclude that pressure tries to resist formation of NS. 
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B. Study of Geodesies 

For simplicity of calculation here we shall consider as before the marginally bound case 
(f( r ) = 0) with n = 3/2. Then R(t,r) has the explicit solution (choosing the initial time 
U = 0) which can be written in convenient form as 



t(r) 



2 

50) 



C(r) - \g{r)B?l* - ^C(r)-^g(r)r^ 



(66) 



To examine whether the singularity at t = to,r = is naked or not, we investigate 
whether there exist one or more outgoing null geodesies which terminate in the past at the 
central singularity. In particular, we shall concentrate to radial null geodesies only. 

First we assume that it is possible to have one or more such geodesies and we choose 
the equation of the outgoing radial null geodesic (ORNG) which passes through the central 
singularity in the past as (near r = 0) 

toRNG = t + a r c (67) 

to leading order in the (t, r) plane with a > 0, £ > 0. 

Now the expression for the singularity time (characterized by R(t s (r),r) = 0) from (66) 



is 



ts(r) = 



y/W)-\jc{r)-\g{r)r^ 



(68) 



and hence the time for central singularity is 



to = 



30 



Here we choose for C(r) and g(r) as 



Co — \ Co 



t3o 



(69) 



C(r) = Car 3 + C k r k+3 
g(r) = g r 3 / 2 + g t r l + 3 / 2 



(70) 



where Co, go & re constants and Ck{< 0) and g t (< 0) arc the first non-vanishing term 
beyond Co and go respectively. As a consequence the expression for t s (r) becomes 



(4- — 

90 \ VCo ^ 



1 



!3o. 



50 \3^fCo~ 
we shall now study the following possibilities: 



§3o 



a/Co — \ Co — §3o 



flo 



r l +.. 



(71) 



(i) k <l , (ii) k > I 



Case I : k < I 

Here for t s (r) we write 

t s (r) =t - 



Ck 

go 



'Co - §3o 



(C fc < 0) 



(72) 
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Now comparing with the geodesic equation (67) we get the relations 



(a) £ > k or (b) £ = k and a < — 



9o 



When £ > k then near r = the solution for R simplifies to 



R = r 



3 2 3 / [Z 2 C fc r- fc 



lsJCo~ 



\9o. 



2/3 



Further for the given metric an ORNG should satisfy 



dt 
dr 



R' + Rv' 



(73) 



(74) 



(75) 



To examine the feasibility of the null geodesic starting from the singularity, we combine 
equations (67) and (74) in equation (75) and we get (upto leading order in r) 



4-1 



2k 



3Cfeto 



l^Co~ 



ho 



2/3 



r~ , (!/_, ^ 0) 



which implies 



2k If 2k 

£ = 1 + y and a = - I 1 + u_ r + — 



3Ckto 



As £ > fc, so from (77) k < 3. Since fc is an integer, we could have 



2/3 



(76) 



(77) 



k = l, £=§ 



or 



On the other hand for £ = k, as before we get fc = 3 and 



(78) 



-1/3 



31, / / 2~ \ (2 + 1/ , +^)C k t 
■- 1(1 + i/_J a 50 t + 2a^C - - 9o + A ^ y 



(79) 
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Case II: k > I 

In this case 



2 fi, /7T 9o - 3C 



t s (r) =t -^-\VCo+ - r ' 



9 ° V Z\fc* Z l9o, 
Now matching the geodesic equation as above we get 



(80) 



(a) £>l or (b) £ = l and a < 

9o 



'Co 



9o ~ 3C 



bo. 



Then for £ > Z we get 



i = 1, £ = I or I = 2, £ = I 



#90 



2/3 



and for £ = I 



(81) 



(82) 



/ = 3, a 



3 3/2 

- (g tl + 2at g ) - - a J C - -go + 



9ito 



\9o 



-1/3 



3 f f 21 

1 + u_ 1 + — ) 50*0 + 2ai (l + v_ 1 )go 



3 M , ./I 2~ , (l + t /_ 1 + f) g| t 

2 v 3 2^6-0 - y 

(83) 



We note that the expressions for 'a' is very complicated both in equations (79) and (83). 
So no definite conclusion is possible on the role of pressure in determining in the final state 
of collapse by ORNG. 



VI. DISCUSSIONS AND CONCLUDING REMARKS 

An extensive analysis of the four dimensional Szckeres model has been done for the matter 
containing pressure. When matter is in the form of perfect fluid then the isotropic pressure 
turns out to be a function of time only while the matter density is a function of all the four 
space-time variables. In this case, assuming a polynomial form for pressure, cosmological 
solutions have been obtained and their asymptotic behaviour have been studied. Both in 
quasi-spherical and quasi-cylindrical model the solution approaches isotropy along fluid 
world line as t — > oo. 

Secondly, for the matter with tangential stress only, solutions are possible for quasi- 
cylindrical model. Here both the tangential stress and the matter density turns out to be 
a function of t and r only. The scale factor R has parametric solution as for dust model 
and does not depend on the tangential stress. However, choosing the parameter C\ =0, R 
has a power law solution and it is possible to have a complete solution if we assume the 
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tangential stress proportional to t 2 . 

Lastly, gravitational collapse has been studied in details for anisotropic pressure (i.e., 
both radial and tangential pressures are non-zero and distinct) in quasi-spherical model. 
Here we have to assume the radial pressure as a function of r and t of the form (see eq. 
(50)) p r — g(r)/R n . Also to solve the differential equation in R (see eq. (53)) we consider 
only the marginally bound case (i.e., / = 0) only. Then equation (64) shows a comparative 
study between the time of formation of trapped surface and the time of formation of 
central singularity. To simplified further, we choose n = 3/2 and detailed analysis has been 
done using equation (65). Table I shows all possibilities for the parameters involved in 
the expression. If the initial density gradient at the centre is positive definite (or negative 
definite) then as in dust case we have definitely a black hole (or naked singularity) as the 
final state of collapse. But when p\ has no definite sign (as in third and fourth cases) then 
for black hole solution it is possible to have negative density gradient at the centre initially. 
In fact, near the singularity if the initial density and pressure has identical behaviour (i.e., 
increase or decrease simultaneously) then even with negative density gradient (initially at 
the centre) we can have black hole as the end state but if the initial density and pressure 
has opposite behaviour (i.e., one decrease while other increases and vice versa) then we 
have identical character as in dust case. So we conclude that pressure tries to resist the 
formation of naked singularity. Finally we have studied the geodesies to examine whether 
it is possible to have any future directed non space-like geodesic terminating in the past at 
the singularity. For simplicity, we have considered only radial null geodesic and it is found 
that the end state of collapse is characterized by the coefficients of the series expansion of 
initial density and pressure (radial). Due to complicated expressions we can not definitely 
characterize the role of pressure. Therefore, in the context of local visibility, we say that 
pressure tries to cover the singularity. 
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